c1=new TCanvas
  gStyle->SetOptStat(0)

  c1->SetTheta(34.76695);
c1->SetPhi(57.67241);
c1->SetFrameBorderMode(0);
  
h3rawmean->SetXTitle("log_{10}#beta#gamma");
h3rawmean->SetYTitle("Timebin");
h3rawmean->SetZTitle("<Q_{cluster}'>_{60%} (a.u.)")
h3rawmean->SetTitle("Timebin Gain Calibrated");



h3rawmean->GetXaxis()->SetTitleOffset(2.2)
h3rawmean->GetYaxis()->SetTitleOffset(1.8)
h3rawmean->GetZaxis()->SetTitleOffset(1.4)
h3rawmean->Draw("surf1")

c1->Print("3dph_tbQcalib.C")
c1->Print("3dph_tbQcalib.png")

//================

    c1=new TCanvas
  gStyle->SetOptStat(0)

  c1->SetTheta(34.76695);
c1->SetPhi(57.67241);
c1->SetFrameBorderMode(0);
  
h3rawrms->SetXTitle("log_{10}#beta#gamma");
h3rawrms->SetYTitle("Timebin");
h3rawrms->SetZTitle("RMS_{60%}(Q_{cluster}') (a.u.)")
h3rawrms->SetTitle("Timebin Gain Calibrated");



h3rawrms->GetXaxis()->SetTitleOffset(2.2)
h3rawrms->GetYaxis()->SetTitleOffset(1.8)
h3rawrms->GetZaxis()->SetTitleOffset(1.4)
h3rawrms->Draw("surf1")

c1->Print("3dph_tbQcalib_rms.C")
c1->Print("3dph_tbQcalib_rms.png")

//==============
c1=new TCanvas
  gStyle->SetOptStat(0)

    
  h3rawmean->ProjectionY("ph_41",41,41)
  ph_41->Draw("hist")
  ph_41->SetYTitle("<Q_{cluster}>_{60%} (a.u.)")
  ph_41->SetTitle("#beta#gamma=10^{3.5}, Timebin Gain Calibrated")
  ph_41->SetXTitle("Timebin")
  ph_41->Draw("hist")
  c1->Print("1dph_tbCalib.png")
  c1->Print("1dph_tbCalib.C")

  //==============
c1=new TCanvas
  gStyle->SetOptStat(0)

  htbratio->SetMinimum(0.4);
htbratio->SetMaximum(0.8);
htbratio->Draw("colz")
htbratio->SetXTitle("log_{10}#beta#gamma");
htbratio->SetYTitle("Timebin");
htbratio->SetTitle("#epsilon (Timebin Gain Calibrated)");
c1->Print("epsilon.png")
c1->Print("epsilon.C")
    
//================

    c1=new TCanvas
  gStyle->SetOptStat(0)

  c1->SetTheta(34.76695);
c1->SetPhi(57.67241);
c1->SetFrameBorderMode(0);
  
h3trumean->SetXTitle("log_{10}#beta#gamma");
h3trumean->SetYTitle("Timebin");
h3trumean->SetZTitle("<Q_{cluster}'> (a.u.)")
h3trumean->SetTitle("Timebin Gain Calibrated");



h3trumean->GetXaxis()->SetTitleOffset(2.2)
h3trumean->GetYaxis()->SetTitleOffset(1.8)
h3trumean->GetZaxis()->SetTitleOffset(1.4)
h3trumean->Draw("surf1")

c1->Print("3dphTru_tbQcalib.C")
c1->Print("3dphTru_tbQcalib.png")

 
//================

    c1=new TCanvas
  gStyle->SetOptStat(0)

  c1->SetTheta(34.76695);
c1->SetPhi(57.67241);
c1->SetFrameBorderMode(0);
  
h3trurms->SetXTitle("log_{10}#beta#gamma");
h3trurms->SetYTitle("Timebin");
h3trurms->SetZTitle("RMS(Q_{cluster}') (a.u.)")
h3trurms->SetTitle("Timebin Gain Calibrated");



h3trurms->GetXaxis()->SetTitleOffset(2.2)
h3trurms->GetYaxis()->SetTitleOffset(1.8)
h3trurms->GetZaxis()->SetTitleOffset(1.4)
h3trurms->Draw("surf1")

c1->Print("3dphTru_tbQcalib_rms.C")
c1->Print("3dphTru_tbQcalib_rms.png")
